Continuous melting of compact polymers 
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The competition between chain entropy and bending rigidity in compact polymers can be ad- 
dressed within a lattice model introduced by P.J. Flory in 1956. It exhibits a transition between an 
entropy dominated disordered phase and an energetically favored crystalline phase. The nature of 
this order-disorder transition has been debated ever since the introduction of the model. Here we 
present exact results for the Flory model in two dimensions relevant for polymers on surfaces, such 
as DNA adsorbed on a lipid bilayer. We predict a continuous melting transition, and compute exact 
values of critical exponents at the transition point. 

Condensed phases of polymers are ubiquitous in the hving world. Examples range from the organi- 
zation of DNA in viruses, to the native state of globular proteins. For example, DNA condensation is 
required in order to fit the complete viral genome, which can be tens of microns long, into the viral capsid 
typically tens of nanometers in diameter. As the persistence length of the DNA (50 nm) is comparable 
to the linear dimension of the viral capsid (10-100 nm), one would expect it to pack into ordered, crys- 
talline configurations. Indeed, cryo-electron-microscopy studies of the T7 virus have revealed a circular, 
inverse-spool arrangement of DNA in the capsid 0| . The observed condensed state of DNA is a result 
of many competing effects such as chain and solvent entropy, electrostatic interactions, bending stiffness 
of the polymer, and interaction with capsid proteins . 

These observations motivate the study of the competition between chain entropy and bending energy 
of a polymer chain under extreme confinement, when only compact configurations are allowed. The 
resulting transition between an entropy dominated disordered state and an energetically favored crys- 
talline one is captured by a simple lattice model of a compact semiflexible polymer, first introduced by 
Flory almost 50 years ago in the context of polymer melting 0] . It has also served as the starting point 
for constructing simple lattice models of protein folding ^] . 

In this Letter, we provide an analytic solution of the Flory model in two dimensions using field theo- 
retical methods. We resolve the long-standing debate [MSHIalMllfl about the nature of the transition, 
where some numerical results and approximate analytic approaches |^ were indicative of a first-order 
transition, while other numerical work ll(l | and theoretical arguments based on similarities with the 
exactly solvable six- vertex model 0, 0| , seemed to imply a continuous transition. We show that the 
melting transition of a single semiflexible compact polymer is continuous. Furthermore, we characterize 
the transition precisely by calculating the exact values of the critical exponents that accompany it. 

Semiflexible loop model In the Flory model a configuration of a polymer chain is described by a 
random walk on a square lattice (cubic, in three dimensions) which visits every site of the lattice exactly 
once. The energy associated with a configuration is Be > 0, where B is the number of bends in the 
chain. At large temperatures {ksT 3> e) bends proliferate for entropic reasons, thus reducing the 
persistence length of the polymer. In the low-temperature limit {ksT <C e) bends are very unlikely and 
the persistence length is large, of the order of the lattice size, as seen in fig. ^ The central question 
we address here is the nature of the transition between the disordered high-temperature phase and the 
ordered low-temperature phase. 

The interplay between the entropy and the energy is encoded in the partition function of the Flory 
model 

= (1) 

where w = exp(— e/ZcsT) is the Boltzmann weight associated with a bend, and BlTp) is the number of 
90° bends in the compact chain configuration F^. For T — > (w — > 0) the equilibrium state is one that 
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FIG. 1: Phase diagram of the Flory modeL The Bohzmann weight of a single bend in the chain is w. For ix; = 1 
the model describes a flexible compact polymer with many bends present in a typical polymer configuration. 
For w < 1 there is an energy penalty associated with bending and we are dealing with a semiflexible chain. For 
Wc < w < 1 the chain is critical. At w = Wc there is a continuous transition to an ordered chain configuration; the 
order parameter is the difference between the number of vertical and horizontal links. The polymer configurations 
shown in the top three panels (with periodic boundary conditions imposed in the horizontal direction) are typical 
for w — 1/3, 10 = Wc ~ 1/2, and w = 1. The density plots below them are of the second height component in 
the related loop model, for the same values of w, clearly showing the transition from a rough to a flat interface. 



minimizes B, while for T oo (w —^ 1) entropy dominates and there is a large, macroscopic number 
of bends in a typical chain configuration; see fig. ^ In order to investigate the transition between the 
high and low temperature phase we construct a mapping of the polymer problem to a lattice model of 
fluctuating loops, following Nienhuis 

We define the semiflexible loop model (SFL), whose configurations Tl are all possible ways of covering 
the square lattice with non- intersecting loops. The partition function of the SFL model reads 

ZL=Y.n''^^-^w^^^-\ (2) 

where each of the N{rL) loops in the configuration Tl carries a weight n > 0, and bends are weighted 
by w as in the Flory model. Note that the n — > limit of / n yields the Flory model partition function 
for a single closed compact polymer. Instead of explicitly computing the approach we take is to 
construct a continuum limit of the SFL model in the form of a field theory. Computations of correlation 
functions in this field theory lead to the phase diagram and exact scaling exponents for any n, including 
the polymer (n — > 0) limit. 

Height model Since loops in F^ do not intersect they can be interpreted as contour lines of the height 
of an interface. In the SFL model, the height is a three component scalar field. To see this, we first note 
that bonds not covered by loops form so-called ghost loops, if periodic boundary conditions are imposed. 
Then we assign to each real and ghost loop an orientation (clockwise or anticlockwise) independently and 
randomly. By doing so, each bond of the square lattice is found in one of four states, labeled A, B, C, 
or D, depending on whether it is covered by a real or a ghost loop oriented in one or the other direction. 
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We parametrize n = 2cos(7re), with < e < 1/2, and assign a weight exp(±i7re) (resp. exp(±i7r/3)) 
to each orientation of a real (resp. ghost) loop. Tracing over the two loop orientations, we recover the 
desired weights: n for the real loops and 1 for the ghost loops. 

The heights for the SFL model are defined at the centers of the elementary square plaquettes of the 
lattice. Height differences between neighboring plaquettes are given by vectors A, B, C, or D. Since each 
vertex is shared by four bonds, all necessarily in different states, a single relation, A + B + C + D = 0, 
results from requiring that the total height increment when traversing a closed lattice path be zero. Thus, 
we conclude that out of the four vectors three are linearly independent and that the heights are three 
dimensional. For the remainder of the paper we adopt the following normalization: A = (— 1,+1,+1), 
B = (+1, +1, -1), C = (-1, -1, -1), and D = (+1, -1, +1). The partition function of the SFL model 
can now be rewritten as a sum over height configurations, since a height configuration uniquely specifies 
an oriented loop configuration. 

The continuum limit of the SFL model is obtained by coarse-graining the height field over domains 
much smaller than the system size and much bigger than a single plaquette. This yields a continuous 
height function h(x) = (ft,^(x), ft,^(x), /i'^(x)) defined over the plane x. For n < 2 the height field is 
assumed rough |l2l | and the partition function can be written as a sum (path integral) over all height 
functions with weights given by exp(— S*), where 

= i ^d^x K^j{dh' ■dh^) + ^ J d^x (eo ■ h)p(x) +7 yd^x cos {2Trh^) . (3) 

This particular form of the coarse-grained free energy S is dictated by the symmetries of the loop model 
which involve translations in the space of heights, and permutations of the height-difference vectors [l^. 
The value of the background charge 

eo = -|(l/3 + e,0,l/3-e) (4) 

follows from imposing the correct weighting of loops that close at infinity (p(x) is the curvature of the 
x-surface), while the non-zero elastic constants 

i^ii = J(5/3-e), ifi3 = ^(e-l/3), ^22 = ^^22(6,1^') (5) 

are partially fixed by the loop ansatz which states that loop weights do not flow under renormalization 

In the previously studied infinite-temperature limit {w —>■ 1), the elastic constant K22 = 7r(l — e)/(5 — 
3e) is also fixed by the loop ansatz, and the cosine term is absent from S -f'or w < I the cosine 

potential accounts for the energy penalty associated with bends. Its scaling dimension in the w = 1 
theory, xb = t^/K22 = (5 — 3e)/(l — e), is greater than two for < e < 1/2 (0 < n < 2), and is therefore 
an irrelevant perturbation. It has the effect of renormalizing the elastic constant 7^22 which ultimately 
drives the transition to the crystalline phase at w ~ Wc- 

Phase diagram and critical exponents The field theory of the SFL model implies the following phase 
diagram for any n < 2, including the polymer limit n — > 0. For w = I the interface described by the 
height field is rough and the loop model is critical, i.e., it is characterized by a power- law distribution of 
loop sizes. As w is reduced, the bending stiffness of the loops increases, and the height remains rough 
until a critical value w = Wc- At there is a roughening transition jl4j to a smooth interface; see 
fig.n In the loop model this corresponds to a phase with a vanishing density of bends. The roughening 
transition is continuous, in the Kosterlitz-Thouless universality class. At the transition the cosine term 
in S becomes relevant in the renormalization group sense jl4l |. In other words, the scaling dimension of 
the operator cos(27rft,^) at w = is 2. Computing this dimension using the field theory defined by S 
then leads to the prediction K22{e,Wc) — n/2. One important test of this result comes from the n — 1 
(e = 1/3) SFL model, which maps to the exactly solvable six- vertex model. In this case = 1/2 and 
^22 = arcsin( 1/210) have been computed previously [l5| . 

In the polymer limit the transition described by the field theory takes us from a critical geometry of 
a disordered compact polymer, at low bending stiffness, to an ordered state, at high bending stiffness; 
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7r/2 


-5/8 


-0.7 



TABLE I: Numerical transfer matrix results. Exact (for w = l,Wc) and numerically determined values of the 
elastic constant 7^22 were used to compute the exponent 6th from eq. for various values of the bend weight 
w; this compares favorably with the direct numerical determination Onum- 



see fig. n The critical phase is characterized by scaling exponents ly and 6. Namely, the probability 
distribution p{r,l), for the end-to-end distance r of a polymer of length has a scaling form p(r, I) = 
r^fir/l") IfiJ- Using the field theory defined by eq. lO these exponents can be computed exactly |l2l| . 
We find 

1 , . 3^ - 16i^22 

1^ = -, and 6= ; (6) 

z on 

note that 9 varies continuously with w via the elastic constant K22- The calculation of i> gives another 
important check on the field theory since v = 1/2 follows directly from the fact that compact polymers 
completely fill the plane. Finally, exact values for K22 in the w ~ 1 and w = Wc case, lead to predictions 
d{w = 1) = 5/56 and d{w = Wc) = —5/8. These values of 9 have a simple physical interpretation. 
Namely, in the absence of bending stiffness {w = 1) the ends of the chain repel (positive 9), as is the 
case for open chain conformations of polymers in a good solvent. Somewhat surprisingly, this entropic 
repulsion decreases and turns to attraction (negative 9) as the chain stiffness is increased. 

Numerical transfer matrix results The main result of this paper, the prediction of continuous melting 
of compact polymers in two-dimensions, rests on the validity of the proposed field theory of the Flory 
model, eq. To further test the field theory we have made use of a numerical transfer matrix approach 
[1^ to compute the elastic constant K22 and the exponent 9. Results are shown in table We see that 
eq. (O is satisfied to a very good approximation and that the exact values for the two quantities at 
w = 1 and w = Wc are confirmed by the numerical computations. 

Experiments DNA absorbed on a lipid bilaycr provides a laboratory for testing theories of two- 
dimensional polymers. Previous experiments 17] on this system have measured the theoretically pre- 
dicted scaling exponents derived from the two-dimensional self-avoiding walk model. These measure- 
ments were done at low DNA concentrations. As the amount of DNA on the surface is increased it was 
observed that single chain conformations become more compact. To study the transition discussed here 
compaction to areas of linear dimensions approaching the DNA persistence length would be needed. It 
would then be interesting to monitor the two ends of the DNA using fluorescent labels, as the amount 
of compaction is increased. This would provide direct information about p(r, /) and the exponent 9 
computed above. 
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